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Abstract 

Aim of this paper is to give an extensive treatment of bipartite mean field spin systems, 
ordered and disordered: at first, bipartite ferromagnets are investigated, achieving an expUcit 
expression for the free energy trough a new minimax variational principle. Furthermore via 
the Hamilton- Jacobi technique the same free energy structure is obtained together with the 
existence of its thermodynamic limit and the minimax principle is connected to a standard 
max one. 

The same is investigated for bipartite spin-glasses: By the Borel-Cantelli lemma a control of 
the high temperature regime is obtained, while via the double stochastic stability technique 
we get also the explicit expression of the free energy at the replica symmetric level, uniquely 
defined by a minimax variational principle again. 

A general results that states that the free energies of these systems are convex linear combi- 
nations of their independent one party model counterparts is achieved too. 
For the sake of completeness we show further that at zero temperature the replica symmetric 
entropy becomes negative and, consequently, such a symmetry must be broken. The treat- 
ment of the fully broken replica symmetry case is deferred to a forthcoming paper. As a 
first step in this direction, we start deriving the linear and quadratic constraints to overlap 
fluctuations. 

Keywords: equilibrium statistical mechanics, bipartite systems, spin glasses 

1 Introduction 

The investigation by statistical mechanics of simple and complex mean field spin systems is ex- 
periencing a huge increasing interest in the last decades. The motivations are at least two-fold: 
from one side, at the rigorous mathematical level, even though several contributions appeared 
along the years (see for instance [18] [23][32j[2]), a full clear picture is still to be achieved (it 
is enough to think at the whole community dealing with ultrametricity in the case of random 
interactions as in glasses |28| |29|[3]): at the applied level, these toy models are starting to be 
used in several different contexts, ranging from quantitative sociology [7]|13)[T2] to theoretical 
biology [30] [6]. 

It is then obvious the need for always stronger and simpler methods to analyze the enormous 
amount of all the possible "variations on theme", the theme being the standard Curie- Weiss 
model (CW) |^|19| for the simpler systems, or the paradigmatic Sherrington-Kirkpatrick model 
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(SK) [28] [ig for the complex ones. 

As a result, inspired by recent attention payed on two groups in interactions (i.e. decision making 
processes in econometrics [25] [14] or metabolic networks in biology [26] [27]), we decided to focus, 
in this paper, on the equilibrium statistical mechanics of two bipartite spin systems, namely the 
bipartite CW and the bipartite SK. 

At first we approach to the problem of bipartite model by studying in Section 2 the bipartite 
ferromagnet, obtaining, both via a standard approach and through our mechanical interpretation 
of the interpolation method |il7j [4j^l5j [8j. a variational principle for the free energy in thermo- 
dynamic limit. 

In Section 3 we open the investigation of the bipartite spin glass model, and following the path 
yet outlined in [10] [9] we get the annealed free energy (more precisely, the pressure), with a 
characterization of the region of the phase diagram where it coincides with the true one in the 
thermodynamic limit, and the replica symmetric free energy, by the double stochastic stability 
technique, which stems from a minimax variational principle, whose properties are discussed too. 
Finally we calculate the zero temperature observable and, by noticing that the entropy is nega- 
tive defined, we conclude that replica symmetry must be broken. Despite a full replica symmetry 
breaking scheme deserves a whole work, here we start introducing its typical linear and quadratic 
constraints, obtained by Landau self-averaging of the internal energy [16]. 
A last section is left for conclusions and outlooks. 

2 Bipartite ferromagnets 

We are interested in considering a set of N Ising spin variables, in which it is precisely defined a 
partition in two subsets of size respectively Ni and N2. We assume the variable's label of the first 
subset as ai, i = 1, iVi, while the spins of the second one are introduced by Tj, j = 1, A^2- 
Of course we have A^i + = N, and we name the relative size of the two subsets Ni/N = a^, 

N2/N = 1 - OAT. 

For the sake of simplicity, in what follows we deal with both parties formed by dichotomic 
variables, while we stress that the method works on a very general class of random variables 
with symmetric probability measure and compact support |15| . 
The spins interact via the Hamiltonian Hj\f-^^_]\f^{a,T,hi,h2), with hi > 0, /12 > 0: 

2 A^i N2 Ni N2 

HnuN2 {(^, T, hi , /12) = -J^—^ 5Z ^^■^i - hl^fTi- h2Y^ Tj. 

^ i=l j=l 1=1 j=l 

We notice that spins in each subsystem interact only with spins in the other one, but not 
among themselves; we have chosen to skip the self interactions to tackle only genuine features 
stemming from the exchange ones. The reader interested in a (different) treatment of bipartite 
ferromagnetic models with self interactions may refer to [14]. 

Partition function Z, pressure A and free energy per site / are defined naturally for the model: 

Z^,^N2W,hi,h2) = ^^e-'3^^i>^2(-.-''^i.'^2)^ 

<T T 

AN^^N2{l3,hi,h2) = ^ logZNi,N2{(3,hi,h2), 

iVi + iV2 

fNi,N2{f3,hi,h2) = ——ANi,N2i(3,hl,h2), 
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while the thermodynamic hmit of A, f will be denoted via A{a, /3, /ii/i2) = limjv A^j^^N^iPj hi, h2) 
and /(a, /3, /ii, /i2) = limTv fNi,N2Wj hi, /12), where wc stressed the prescription adopted in taking 
the infinite volume limit, performed in such a way that when N, Ni, N2 — s- 00, Ni/N — t- a € (0, 1), 
and N2/N ^ 1 — a G (0, 1). Taken z{a,T) as a generic function of the spin variables, we can 
also specify the Boltzmann state of our system as 

^^(^ ^ Eg Et ^(^' ^) exp(-/3gjv,,jv, (0-, T,hi,h2)) 

Zni,N2{P, hi,h2) 

As usual the order parameters (the respective magnetizations of the two systems) are 

1 

i 

N2 

N2 



thus the Hamiltonan reads off as 



Hni,N2 {<^, r, hi , h2) = -N 



2-—^ — mNinN2 + hiuiNi + h2aNnN2 

1 + OiN 



(4) 



Remark 1. As will he clear soon, the choice of the factor 2 in the Hamiltonian is made in such 
a way that the balanced bipartite model with a = 1/2 has the same critical point of the one party 
model, i.e. j3 = 1. 

2.1 The occurrence of a minimcLX principle for the free energy 

Now we give the explicit form of the pressure of the model, together with some interesting 
properties. The main result is the following: 

Theorem 1. In the thermodynamic limit, the pressure of the bipartite ferromagnetic model is 
given by the following variational principle 

A{a,P,hi,h2) = m.eo^TCLiTLAtriai{fh,n), (5) 

m n 

where 

■^trial{'rn, n) = log 2+a log cosh (2/3(1 — a)n + hi) + {\—a) log cosh {2Pafh + h2) — 2j3a{l—a)rhn. 

Furthermore, the solution is uniquely defined by the intersection of 

m = tanh (2/3(1 — a)n + /ii) , (6) 
n = tanh (2/3am + /12) , (7) 

for 771 > and n> 0. 
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Proof. The proof can be achieved in several ways (i.e. a direct approach by marginahzing the 
free energy with respect to both the parties); we chose to foUow the path outhned in [9] as this 
may act as a guide later, dealing with frustrated interactions. 

For the sake of convenience, we rename /3hi — )• hi and /?/i2 — )• /12 as switching back to the original 
variables is straightforward in every moment but this lightens the notation. 
To our task we need to introduce two trial parameters m, n mimic the magnetizations inside 
each party, one interpolating parameter t G [0, 1] and an interpolating function A{t) for the free 
energy as follows 

A{t) = - log 5] ^ e* ('''^-(-)" W) (2/^(1--)^ -'+2^°'^ ^T' -0 ^h, ^.+h, E, r, ^ 

U T 

(8) 

such that, for t = 1 our interpolating function reduces to the free energy of the model (the 
pressure strictly speaking), while for t = reduces to a sum of one-body models whose solution 
is straightforward. 

We can then apply the fundamental theorem of calculus to get the following sum rule 

^(1) = A(0) + jjl' ^dt. (9) 

To quantify the latter we need to sort out both the streaming of the interpolating function as 
well as its value at t = 0, namely 

dtA{t) = 2(3a{l — a)[{mn) — n{m) — m{n) + mn] — I3a{l — a)mn, (10) 
A{t = 0) = log2 + alogcosh[2/3(l - Q)n + /ii]+)(l - a)logcosh[2/3am + /12], (11) 

where in equation (jlOp we added and subtracted the term /3q:(1 — a)fhn so to write explicitly 
the sum rule in terms of a trial function Atriai{oi, (3, hi,h2) and an error term 5(m, n): 

A{a,(3,hi,h2) = Atriai{a, P,hi,h2) + S{fh,n), (12) 

where 

Atrial {oi, P, hi, = log 2 + a log cosh[2/3(l — a)n + hi] + (1 — a) log cosh[2/3am + /i2] 

- 2/3a(l - a)mn, (13) 

S{m,n) = /3a(l — q) / {{m — fh){n — n))tdt, (14) 

Jo 

for every trial functions m, n. 

Note that the averages {.)t take into account that the Boltzmannfaktor is no longer the standard 
one introduced in eq.(l), but incorporates the interpolating structure tuned by the parameter t. 
We stress that, at this stage, the error term S{m, n) is the source of the fluctuations of the 
order parameters (which are expected to reduce to zero in the thermodynamic limit and give the 
label S) is not trivially positive defined as in many other cases, even hardly to investigate (i.e. 
mono-party spin-glasses [ISj), however the idea of choosing properly m, n so to make it smaller 
and smaller (eventually zero) still holds obviously. 

So we study at fixed n the behavior of our trial function in fh by looking at its derivative: 

drhAtriai{c(,(3,hi,h2) = 2/3a(l - Q)[tanh(2/3am -1- hi) - n], (15) 
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so, at given n it is increasing in m and AtriaiiP, /ii, ^2) convex in m. 
By a direct calculation we get the same result inverting n <^ fh: 

dfiAtriai{a,P, hiM) = 2/3a(l - a)[tanh(2/3(l - a)n + /12) - fh]- (16) 

As it is crystal clear that the roles of m, n are of the local magnetizations, we can allow ourselves 
in considering only values fh > tanh(^/ii) such that there exist a unique n(m) > : tanh[2/3(l — 

a)h{fh) + hi] = fh. 

From now on let us switch from Atriai{(^, P, hi, /12) to a less rigorous labeling Atriaii^n, h) which 
aims to stress the relevant dependence by its variables time by time: For this value h{fh) lastly 
obtained, the trial function has its minimum in n at fixed fh such that we can substitute it and 

get Atrial{fh) = Atrial{fh,h{fh)). 

Now, as 

dmAtriai{fh) = 2/3a(l — a)[tanh(2/3am — hi) — h{fh)], (17) 

we can consider A[fh) as a function of tti^ to see easily that A is concave in tti^, so it has its 
unique maximum where its derivative vanishes. 

Overall we can state that A{a, hi , ^12) = sup^ inf j^ A{fh, h), whose stationary point is uniquely 
defined by the solution of the system of self-consistence relations 

fh = tanh (2/3(1 - a)n + /ii) , (18) 
n = tanh {l^afh + /i2) , (19) 

as stated in Theorem 1. 

Furthermore, if we restrict ourselves in considering h = h{fh) -as imposed by the variational 
principle- the error term (the fluctuation source) results positive deflned: This statement can be 
understood by marginalizing with respect to the r party the free energy (summing over all the 
T-configurations) so to substitute n(r) by tanh[2/3am(cj)t + 2/3afh{l — t) + hi] and noting that 
m(cr) > fh implies tanh[2^ami -|- 2j3afh{l —t) + hi] > h{fh) such that 

{{m — fh)[tanh.{2a/3fh + h2) — n] > 0, 

the error term is positive defined. 

Now, in order to show that the error term is zero in the thermodynamic limit (such that the 
expression of the trial becomes correct) we proceed on a different way: the idea is to marginalize 
with respect to one party, so to remain with a single ferromagnetic party with a modified external 
field and then use the standard package of knowledge developed for this case. 
So, at first, we marginalize the free energy with respect to the r variables: 

ANi,N2i/3,hi,h2) = ] logZN,,N2{(^,hi,h2) = (20) 

iVi + 1\2 

= ^^^^^log5^2^2cosh^2(2a^m + /i2)exp(/n^(7,)= (21) 

= (1 - a) log 2 + ^- log y log cosh(2/3am+h2)+fe2 Ei . (22) 

^ ' NI + N2 



a 



We can use now the convexity of the logarithm of the hyperbolic cosine as log cosh (re) > 
log cosh(x) + {x — x) tanh(x) to get 

g-iV2iogcosh(2a^m+ft2)g^2iogcosh(2a^m+h2) + N22aP{m - m) tanh(2a^m + h2) < 1, 
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by which we bound the free energy through a new trial function Atrial as 

AN,,N,(.(3,h,,h2) = (l-«)log2+ [ iQgy-giV.logcosh(2^am+fe.)+/..E,^i > (23) 

iVi + I\2 ^ — ' 

> {I — a) log 2 + log cxp(iV2 log cosh(2Q;^m + /t2)) x 

JSi + A/2 

X exp(-2a/3A/2"^ tanh(2a^rn + ^2)) exp{2a/3 tanh.{2aPm + /i2) ^ CTj) =(24) 

j 

= (1 - a)(log2 + logcosh(2a^m + /12)) - 2a(l - a)/3mtanh(2Q;;3m + /12) + 
+ a log 2 + a log cosh[2/3(l — a) tanh(2a^m + /12)] = ^tr-iai- 
Once defined h{m) = tanh(2a/3m + /i2), we can look for the fh streaming of the trial, namely 

^rnAtr^al{m) = {2(3af {1 -«)(!- n^) [tanh (2/3(1 - a)h) - fh]. (25) 

If we now consider the streaming with respect to h of Atrial we get 

dhAtriai{m) = {dnfh)drfiAtrial{m) = 2/3a(l - a)[tanh(2^(l - a)n) - m(n)]. (26) 

So the streaming is decreasing in h and the trial is concave in n, there exist a unique maximum 
where the derivative vanishes: Properly choosing m — >■ m and n — >■ n we get the statement of 
the theorem. □ 

As we are going to deepen the extremization procedure in these bipartite models, we stress 
that an important feature that seems to arise from our study, is the occurrence of a minmax 
principle for the free energy, usually given by a maximum principle in the ordered models, and 
a minimum principle in the frustrated ones. 

We finally report an interesting result about the form of the pressure (or equivalently, the free 

energy) of the model. Indeed it turns out to be written as the convex combination of the 
pressures of two different monopartite CW models, at different inverse temperatures, as stated 
by the following 

Proposition 1. In the thermodynamic limit, the following decomposition of bipartite ferromag- 
netic model free energies into convex sums of monoparty ones is allowed: 

A{a, P, hiM) = aA^'^iP', hi) + (1 - a)A^'^{P", ^2), 

with P' = 2^(1 - a) I and P" = 2Pa^. 

Proof. We start setting the trial values for the inverse temperatures of the two monopartite 
models, as /3' = 2/3(1 — 0)0^ and /3" = 2j5a/a^^ with a a real parameter to be determined later, 
and we set the external fields to zero for the sake of clearness. It is 



]mn] 



and since 2mn < rn^o? + Va / 0, we have 

(22 2 
2^iVa(l - a)^ + 2pNa{\ - 
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thus we conclude, when N ^ oo 

A{a,^)< aA^^ (/?') + (1 - a)A^^ (/3") (27) 
The inverse bound is proven noticing that 

Zn,,nM _ exp (2/3Na{l - a)mn - ^ - ^) exp + 

Z^,(/3')Z7V,(/3") " Z,v,(/3')^JV,(/3") 



/3"n2 



> exp $7 



exp ^2/3A^q(1 — a)mn ^ ^ 
2/3iVa(l - a)mn ^ ^ 



exp ( 2/3A^a(l — a)mn 



2 2 
2~ 

where we have denoted with Q the joint state of the two monopartite systems. Now, bearing in 
mind the definition of /?' and f3" given at the beginning, we get 

ZNi^NjjP) ^ -2NI3a{l-a)[am-^Y 

z^.mzN.m - 

and then, in thermodynamic limit 

A{a,l3)>aA^^{l3') + {l-a)A^^{0')-2l3a{l-a){[afh-'^ ) (28) 

Now we notice that the extrema in (|27p and (|28p . respectively a minimum and a maximum, are 
obtained with the choice = ^, that completes the proof. □ 

It is worthwhile to remark that the two monoparty models here are completely independent, 
in the sense that the order parameters are given by the relations 

fh = tanh(/3'm), 
n = tanh(/3"n). 

2.2 The free energy again: a maximum principle 

Our aim is now to front the mathematical study of this model with the approach described in 
|17|[4][T5][8]. based on a mechanical interpretation of the interpolation method. The problem of 
finding the free energy in the thermodynamic limit is here translated in solving an Hamilton- 
Jacobi equation with certain suitable boundary conditions, and an associated Burgers transport 
equation for the order parameter of the model. In order to reproduce this scheme, with the 
freedom of interpretation of the label t for the time and x for the space, let us introduce now the 
(x, i)-dependent interpolating partition function 



ZN{x,t) = / ^/ ^ exp tajv(l — aj^)mi^nN + 

a T 

, (/3 - / 2 2 2 , A - 2 N , ( I -UN . 

' (a a rrij^ + I n^j + xyaaNrriM un) 



+ hiaNiriN + h2{l - aN)nN 
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such that the thermodynamical partition function of the model is recovered when t = 2(3 and 
X = 0. At this level a is a free parameter to be determined later. We can go further and explicitly 
define the function ^ 

ipN{x, t) = log Zn{x, t), (29) 

that therefore is just the pressure of the model for a suitable choice of {x,t). Now, computing 
derivatives of (ffyr{x,t), we notice that, putting = aaNrri]^ — njv, it is 

dt(pN{x,t) = -^{D%) {x,t), 
dx(pN{x,t) = {DN){x,t), 
dlipN{x,t) = ^(^{D%)-{D^f 

Thus we can build our differential problems trough an Hamilton- Jacobi equation for ipN{x,t) 

( dtipN{x,t) + \{d^ipN{x,t)f + i^dlipN{x,t) = Q inMx(0,+oo) 
\ (^7v(x,0) = aiv^^f (/?', a(/ii + x)) + (1 - a^)^^f (/3", (1 - a)(/i2 - x)) on R x {t = 0}, 

(30) 

where is the pressure of the Curie- Weiss model made of by A^i a spins with inverse tem- 

perature and A^^ is the same referred to the N2 r spin with inverse temperature , with 
/3' = 2/3o^(l — a) and /3" = 2/3a~^a, and trough a Burgers equation for the velocity field DN{x,t) 

dtDN{x,t) + DN{x,t)d^DN{x,t) + ^d^DN{x,t) = inRx(0,+oo) 
Dn{x, 0) = am{(3, hi + x) — {1 — aN)n{aN(3, h2 — x) on R x {t = 0}, 

This is true of course for every choice of the parameter o, that has the role of balancing the 
weights of the single party contributions. Since we have seen that the function (^at is decreasing 
in time, if we put x = 0, with no external fields, we gain 

An„nM < a^^^r(/3') + (1 - aN)A%Y{(3"). 



If we take = y^-^, such that it is f3' = f3" = /3 = 2fi\J a(\ — a), we have, in the infinite 
volume limit 

A{^) < aA^^i^) + (1 - a)A^^{^), 

that easily give us the critical line of the bipartite model, 2(3y^ a{l — a) = 1, obtained straightly 
by the critical point of the two Curie- Weiss models, (3 = 1. Anyway, for reasons that will become 
clear soon, hereafter it will adopted the different value a = 1. 

Remark 2. We stress that the boundary condition in equation \30(l is always an upper bound for 
ifN- 

In order to work out an explicit solution for the thermodynamic limit of the pressure, in primis 
we notice that the main difference with respect to the single party (namely the Curie Weiss [15j) 
is the more delicate form of the boundary conditions. In fact we have that interactions do not 
factorize trivially (in a way independent by the size of the system). It is 

V9^(x, 0) = a^vA^f (/?', hi+x) + {l- aN)A%^{(3", h2 - x), (32) 



8 



however, it is known [15] how to get a perfect control on the function on the r.h.s. of (f32|) . and 
we have 

^n{x, 0) = a^^^(/3', hi+x) + {l- q)^^^(/3", h2 - x) + O (^^^ . (33) 
Evenly we have for the velocity field in t = 

OiNiT^Niil^' , hi+x) — {l — aN)n]\i2 iP" ^ h2 — x) = am{f3' , hi+x) — {l — a)n{(3" , h2 — x) + ( — 



N 

We obtained an Hamilton- Jacobi equation for the free energy with a vanishing dissipative term 
in the thermodynamic limit, while the velocity field, that is the order parameter, satisfies a 
Burgers' equation with a mollifier dissipative term. 

We stress that our method introduces by itself the correct order parameter, without imposing it 
by hands. 

Remark 3. As the next condition on (p]\}(x,t) is tacitely required by the following lemma, we 
stress that the function Dj\f{x, t) is bounded uniformly in N, a, /?, hi, /12, that implies the function 
ip]sr{x,t) to be Lipschitz continuous. 

We can replace the sequence of differential problems with boundary conditions dependent by 
N with the same sequence of equations but with obvious fixed boundary conditions, that is the 
well defined limiting value of ipN and in t = 0. To this purpose it is useful the following 

Lemma 1. The two differential problems 

r dt(pN{x,t) + l{dcc(pN{x,t)f + ^d'^ipNix,t) = ^nMx(0, +00) 

\ ^Nix,0) = aNA%Yil3',hi + x) + il-aN)A%YW",h2-x) = hNix) onRx{t = 0}, 

(34) 

and 

dt(pN{x,t) + \{dx<:pN{x,t)f + i^dl^N{x,t) = in M X (0,+oo) 

ifN{x, 0) = a^^^(/3', hi + x) + \l- a)A^^(/3", /12 - x) = h{x) on M x {t = 0}, 

(35) 

are completely equivalent, i.e. in thermodynamic limit they have the same solution, ipj\f ^ ip and 
(pN ^ '■P CLfid it is 

\V^N -'fN\<0 

Proof. By a Cole-Hopf transform, we can easily write the general form of 6i\i{x,t) = \ipj\f{x,t) 
If Nix, t)\ as 



log 



/+~dyA(y,(x,t))e-^«-(^) 



f^^dyA{y,{x,t)) 



where we introduced the modified heat kernel A{y, (x, t)) = ^ ^ exp (^—N [(x — y)'^/2t + /i(?/)] ) , 
and iijv(?/) = \h{y) — h]sf{y)\, with lim^v NR^ < 00, Vy. Now we notice that as it certainly exists 
an y* such that 

supi?Ar(y) = y* and lira NR]\f{y*) < 00, 
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hence it is 

SN{x,t) < 1| loge-^^-(^*)| = 1 [NR^{y*)] < ^ (^) , (36) 

that completes the proof. □ 

Of course a similar result holds also for the Burgers' equation for the velocity field D]\f. 
Now the path to follow is clear: the problem of existence and uniqueness of the thermodynamic 
limit is translated here into the convergence of the viscous mechanical problem to the free one. 
We can readapt a theorem that resumes a certain amount of results due to Doughs, Hopf, Lax and 
Oleinik |24||15| which assures the existence and uniqueness of the solution to the free problem: 

Theorem 2. The pressure of the generalized bipartite ferromagnet, in the thermodynamic limit, 
exists, is unique and is given by: 

A{f3, a, hi, /i2) = — 2a(l — a)(3nrh + a log cosh {hi + 2(1 — a)/3n) + (1 — q) log cosh (/i2 + 2af3m) , 

(37) 

where, given the well defined magnetization for the generalized CW model respectively for a and 
T, m{(3,h) and n{(3,h), it is 

m{P,a,hi,h2) = m{2/3,hi-D) (38) 
fi{f3,a,hi,h2) = n{2/3,h2 + D). (39) 

Furthermore it is 

\AN{P,hi,h2) - A{p,a,hi,h2)\ < O (J^^ . (40) 

Proof. Well known results about CW model (see for instance [3]|15)) give us the existence and 
the form of the free solution. We know [15j that the free Burger's equation can be solved along 
the characteristics 

' ' = ' (41) 



X = xo + sD{xq,0), 
where 

D{xo, 0) = Q!m(2/3(1 — a), a{hi + xq)) + (1 — a)n(2/3a;, (1 — a)(/i2 — xq)), 

and it is 

D{x,t) = D{xo{x,t),0) (42) 
= am{2(3{l - a),a{hi +x- tD{xo, 0))) + (1 - a)n{2(3a, (1 - a)(/i2 - x + tD{xo, 0))). 

Then we can notice that 

m{P',hi-x + tD{xo,0)) = tanh {hi - X + t{l - a)n) , (43) 
n{^",h2 + x -tD{xo,0)) = tanh {h2 + X + tarn) , (44) 

which coincide with (I38p and (l39l) when x = and t = 2/3. At this point we know fl5] that the 
minimum is taken for y = x — tD{x, t), such that we have 

f{x, t)(^=o,t=2/3) = [l^^x, t) - ^a^m\^', hi+x- tD{xo,0)) - ^(1 - afn\^", h2 - x 
+ tD{xo, 0)) + Q log cosh {hi-x + m(2/3(l - a) + ta) - t{l - a)n) 



+ {1 — a) log cosh (/i2 + X + tarn — n((l — a)t — 2/3a)) 
= A{(3,a,hi,h2), 



(x=0,i=2/3) 
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where A{(3,a,hi,h2) is just given by (j37p . bearing in mind the right definition of M and N 
Now we must only prove the convergence of the true solution to the free one. But equation (|40[ 
follows by standard techniques, because of the uniform concavity of 

[x-yY , ^ .cw r a u , „a , i-t „,\ aCW , 



2t 



+ aA'-^i(3, hi+y) + {l- a)^^^ (/?, /12 - y) 



with respect to y. In fact we have that, by a Cole-Hopf transform, the unique bounded solution 
of the viscous problem is 

^^{x, t) = 1 log exp [-AT (^^^ + aA^'^il^, h, + y) + {1 - a)^^^(/3, h, - y)) 

and we have, by standard estimates of a Gaussian integral, that 

\ip{x,t) - ipN{x,t)\ <0(^-^ 

i.e. also eq. (140 p is proven. □ 

It is interesting to notice that here the minimax principle discussed in the previous section 
has become a pure maximum principle for the free energy, because of the natural choice of the 
order parameter D, that is in our formalism the analogue of the velocity field. Thus we have 
outlined the framework for stating the next 

Proposition 2. As an alternative to Theorem 1, the free energy of the bipartite ferromagnet can 
he obtained even within a classical extremization procedure as it is uniquely given by the following 
variational principle 

/(a,/3,/ii,/i2) = max - D"^ + a^m'^{D) + {I - afn{Df 

a 1 — 01 
— — log cosh ihi + 13(1 — a)n(D) — log cosh (h2 + Pam(D)) 

f3 13 



3 Bipartite spin glasses 

Let us consider a set of A'^i i.i.d. random spin variables cJi, i = 1, A'^i and let us consider also 
another set of i.i.d. random spin variable rj, j = 1, ■..,N2 = N — Ni. We will consider for sake 
of simplicity only dichotomic spin variables, although our scheme is easily extensible to other 
spin distributions, symmetric and with compact support. Therefore we have two distinct sets 
(or parties hereafter) of different spin variables, and we let them interact through the following 
Hamiltonian: 

i=i j=i 

where the are also i.i.d. r.v., with E[^] = and E[^^] = 1, i.e. the quenched noise ruling the 
mutual interaction between parties. In particular we deal with a J\f{0, 1) quenched disorder. It 
is then defined a bipartite model of spin glass where emphasis is given on its bipartite nature by 
neglecting self-interactions, mirroring the strategy outlined in the first part of the work. 
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For the sake of simplicity, as external fields in complex systems must be considered much more 
carefully with respect to the simple counterparts, we are going to work out the theory neglecting 
them at this stage. 

Of course, once the Hamiltonian is given, it results defined the partition function, the pressure 
and the free energy of the model, as 

^Afi,iV2(/3) = '^exjp{-PHN^^N2{(^,T)) , (46) 

cr,r 

AnunM = j;^Ejlog Zn^nM, (47) 
Jn^nM = -^An„nM- (48) 
We can define also the Boltzmann state for an observable function of the spin variables z{a, r): 
(^Ni,N2{z) = Z^l^^{(3) [2;(a,r)exp {-I3Hni,N2{(^,t))] ■ 

and, as in glasses we need to introduce replicas (equivalent copies of the system with the same 
identical quenched disorder), we can define even the Boltzmann product state as i77Vi,Af2 = 
'-^N-i,N2 ^ ••• ^ '^Ni,N2i where the amount of replicas can be specified time by time. 

Remark 4. In analogy with the prescription introduced in the normalization of the bipartite 
ferromagnet, the factor \/2 in the Hamiltonian is put ad hoc in order to obtain for the balanced 
bipartite spin glass (a = 1/2) the same critical point of the Sherrington- Kirkpatrick one party 
model, as will be clarified in the next section. 

The main achievement of the theory would be a complete control of the free energy, or the 
pressure, in the thermodynamic limit, i.e. for Ni,N2 — )• oo. We stress that two different cases 
may arise: the first is that the size of one party grows faster than the other; the second is that 
the two sizes grows in the same way, such that is well defined the ratio Ni/N — )• a G (0, 1) 
and N2/N — 7- (1 — a) G (0, 1) again for coherence with the strategy outlined in the first part of 
the work and for a general higher interest in this case. We will adopt this latter definition of 
thermodynamic limit, and thus the thermodynamic functions depend by the additional parameter 
a ruling the relative ratio among the parties: 

lim^7Vi,iV2(/3) =^(a,/3). 

We must stress that at the moment no rigorous proof of the existence of such a limit is known 
and that there is a deep connection among this limit and the one of the Hopfield model for neural 
networks [lOj^. Finally we must introduce overlap, that is correlation functions among replicas. 
It naturally arises how in this model we have two kind of such a quantity, one referred to each 
party; in fact we define immediately 



Pab 



N2 

3 
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3.1 High temperature behavior 

We start our study of the model by characterizing the high temperature regime. It turns out 
that the system behaves like the annealed one in a wide region of the (q, /3)-plane, as stated by 
the following 

Theorem 3. The pressure of the bipartite spin glass model does coincide with its annealed one 

AAia,(3) =log2 + f3'^a{l - a) (49) 



in the region of the (a,/?) plane defined by 2/3^-y/a(l — q) < 1. 
Proof. At first we calculate the annealing: 

E[Zn„nM] = exp(iV) (log2 + fail - a)) , (50) 

thus 

AA{a,(3)= lim (N)-' logE[ZN„NM] = log2 + p^a{l - a). 
Ni,N2 

Now, following a standard method |31||10|. we want to use the Borel-Cantelh lemma on ZAr^^jv2/[^A'] 
We evaluate the second moment of the partition function: 



W,t} \ ij 

= ^ exp(^/32^^(2 + (;i2Pi2)) 



^ g2(7V){iog2+/32„(i-a)) ((Ar)2/32a(l - 0)912^12) , (51) 

where we neglected terms leading to an error for the pressure 0{N~'^), replacing a^v with a. 
Now we perform the transformation a — >• aa' and r — )■ rr' in order to get qi2 — )• m and pi2 — ?• n. 
Thus we have 

]gr_^2 1 
E[Zn,,N2Y 

= ^ exp((iVi + iV2)/3' a(l — a)mn) , 

(52) 

with /3' = /3^. The last term, as we have seen in the previous section about bipartite ferromagnetic 
models, is bounded for 2/3' y/aJY—a) < 1, i.e. 2/3^y^a(l — q) < 1, that completes the proof. □ 

We notice that this is a result slightly different with respect to the Hopfield model [lOj. In fact 
we have that the annealed free energy and the true one are the same still at small temperatures, 
depending on the different weights assumed by the parties (that are of course ruled by a). This 
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is reflected by the symmetry of the high temperature region with respect to the hne a = 1/2. 
The following argument is going to clarify this point: 

Given respectively a A'^i x A^i and a N2 x N2 random matrix J and J', with both Jij and 
Jj'j normal distributed random variables for every we introduce the interpolating partition 
function 



a r2t^Ni.N2 c , 2a /2(l-t) v^JVi . (3 2{l-t) v^iV2 r/ 



where o is a parameter to be determined a posteriori. Putting /3' = /3a^\/2a and (3" = 
/3a~^y^2(l — a), and neglecting terms vanishing when N grows to infinity, we can rewrite the 
latter expression as 



Now we introduce the function 



ct>N,,N2 it, /3, /3', n = log Zn,,N2 (/3, /3', t) + J + - a) + A(3^a{l - a)) 



It is easily verified that 

(53) 



M,7V2(t = 0) = ayl^f(/3') + (l-aM^f(/3") 
</.jVi,7V2(t = 1) = ^jVi,7V2(/3) + i(/3"a + /3"2(l-a)-4/32a(l-a)) 

where, of course, A^^ is the pressure for the Sherrington-Kirkpatrick model, that is, the only 
one party model. Furthermore we can take the derivative in t and obtain 



since 2/3'/3"y^ a(l — a) = 4/3a(l — a). Hence we get the following bound for the pressure 

An.M^) > a^^f (/3') + (1 - a)^^f - ^ (/3''a + /3"2(1 - a) + A^^a{l - a)) . (54) 



Now we can fix in such a way that /?' = /3" = /3. As a consequence, it results = {1 — a)/a 
and = 2/3^ Y^a(l — a), and the formula (|54p becomes 

^iVi,7V2(/3) > «^^f (^) + (1 - a)^^f (^) - ^ + /32a(l - a). (55) 

Thus the pressure is always greater than the convex sum of the single party Sherrington- 
Kirkpatrick pressure. The extra term is build in such a way we get an equality in the annealed 
region. In fact we have, if ^5 < 1, i.e. 2/3'^^ya{l — a) < 1, that in thermodynamic limit both 
Af/^ and Af^^ are log 2 + j^'^/A, and therefore we get 

log2 + /32a(l-a) >A{a,f3) >log2 + l3'^a{l-a), 
where, as usual, the upper bound is given by the Jensen inequality. 
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3.2 Replica symmetric free energy 

In order to obtain an explicit expression for the free energy density (or equivalently the pressure 
A{a,f3)), we apply the double stochastic stability technique recently developed in [9]. In a nut- 
shell the idea is to perturb stochastically both the parties via random perturbations; these are 
coupled with scalar parameters to be set a fortiori in order to get the desired level of approxima- 
tion. With these perturbations the calculations can be reduced to a sum of one body problems 
via a suitable sum rule for the free energy; by the latter, the replica symmetric approximation 
can be obtained straightforwardly by neglecting the fluctuations of the order parameters. 
Concretely we introduce the following interpolating partition function, for t G [0, 1] 

ZNt,N2it) = ^^exp{Vt-y= (,ijaiTj) ■ (56) 

cr T ij 

Ni N2 



exp(^/r^[/3V2(l - a)pY,^i^i + PV^^Vjrj]), 



« -J 

where r],f} are stochastic perturbations, namely i.i.d. random variables AA[0, 1], whose averages 
are still encoded into E, and, so far, g, p are Lagrange multipliers to be determined later. 
Now we introduce the interpolating function 

AN,,N2{t) = ^^logZN,,N2{t) + (1 - t)a{l - a)P\l - q){l-p)- 

It is easily seen that at t = 1 we recover the original pressure j4(a,/3), while for t = we obtain 
a factorized one-body problem: 



limAT AAri,Ar2(t = 1) = Aat^^aTj (/3), 

IimjvA7Vi,iV2(i = 0) = Ao(a,/3) + a(l-Q)/32(l-g)(l-p) 



2n .^ (57) 



with 



Ao(a,/3) = -^E log exp(/3 V2(l - a)pY,7]m) 

+ i-Elog^exp(/3y2^^r/,r,) (58) 



= In 2 + aEg log cos ^g^3^/2{l - a)p) + (1 - a) Eg log cosh (5/3^2^), 

where Eg indicates the expectation with respect to the M{0, 1) r.v. g. 

Now we must evaluate the t-derivative of AN-^^N2{t) in order to get a sum rule, namely 

1 



A{t = 1) = Ao{a, /3) + dt [j^A{t)j . (59) 

Denoting via ()t the extended Boltzmann measure encoded into the structure (|56p -that reduces 
to the standard one for t = 1 as it should-, we get three terms by deriving the four contributions 
into the extended Maxwell-Boltzmann exponential, that we call A,B,C and follow: 
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A = ^E-^Y.^^M^^rj)=a{l-a)(3'{l-{quPi2)t), (60) 



N 2VtN 



C = ^E|^j:%.(r,) = -5li^/,=J{l-ta.),). (62) 

i 

So we can build the t-streaming of the interpolant Ajs[^^m2{'^) 

j^An^nA^) = A + B + C-a{l-a)p\l-q){l-p) 

= a{l- a)/32 (1 - (gi2Pi2)t ) - "^^""V ^P (1 - ('Zi2>t ) + 

- - (Pi2)i) - a(l - a)/32(l - q){l - p) 

= -{{q-Qi2){p-Pi2))f (63) 

Now we stress that the Lagrange multiphers q and p can be understood here as trial values 
for the order parameters. According to this point of view, the replica symmetric condition, i.e. 
the request that the overlaps do not fluctuate, is equivalent to impose {{q — qi2){p — Pi2))t — 0) 
bringing us to conclude that in RS regime ^jVi,Af2(0 ^ steady function of t, and then 

yliVi,iV2(i = 1) = ^iVi,iV2(/3) = ^7Vi,iV2(t = 0) = ^o(a,/3) + a(l - q)/32(1 - g)(l -p), 
that is, in the thermodynamic limit 



A{p,q,a,l3) = In 2 + oEg log cosh [gfiyj 2(1 - a)p) + 

+ (1 — ay&g log cos h(5/372ag) + a(l - a)p^{l -q){l-p). (64) 

Now we follow the same considerations exploited in |9j. Indeed, the last expression holds 
barely for every possible choice of the trial values q and p of the order parameters (q'12) and 
{Pi2)- Our purpose is then to fix the right value of q and p, imposed by the RS condition 
{{q ~ Qi2){P — P12)) = 0. In primis we notice that the trial function A, as a function of the trial 
order parameters q, p, is uniformly concave with respect to p. In fact it is easily seen that 

dpA{q,p, a, 13) = a{l - a)/?^ {q - Eg tanh^ {gpyj2{l-a)p 



and since E<;tanh2 (gl3yj2{l - a)pj is increasing in p we have the assertion. Furthermore, for 

any fixed q, the function A takes its maximum value where the derivative vanishes, that defines 
implicitly a special value for p{q): 



q = Eg tanh^ (^g^^2{l - a)p(^)) • 



Of course we have that p is an increasing function of q, with p{0) = 0. Now we are concerned 
about A at a fixed level set, i.e. A{p{q),q), and state that it is convex in q. In fact it is easily seen 
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from the last formula that p{q) jq is an increasing function of g, but of course Eg tanh^ (g/3v^) 
is strictly decreasing, thus 

dqA{q,v) = a{l - a)P^ {p{q) - E^tanh^ {gP,/2^y^ 

is increasing in q, that is equivalent to convexity in such a variable. The last equation specifies the 
right value of the trial order parameter q in the RS approximation. Thus the replica symmetric 
approximation of the pressure results uniquely defined by the minimax principle: 

Theorem 4. The replica symmetric free energy of the bipartite spin glass model is uniquely 
defined by the following variational principle: 

A^''^ {p,q,a, f3) = mmmaKA{p,q,a, /3), (65) 

q p 

where 

A{p, q, a,f3) = In 2 + oEg log cosh {g|3^/2{l - a)p) + 

+ (l-Q)Eglogcosh(5/372^)+a(l-a)/32(l-g)(l-p), (66) 

whose the saddle point is reached at the intersection of the following two curves in the (a, /3) 
plane 

q{a,l3) = Eg tanh^ {g 13 ^/ 2(1 - a)p{a,/3)) (67) 

pia,P) = Eg tanh'^{gf3y/2aq{a,/3)). (68) 

We can go further and put some constraints, imposing the two curves to intersect away 
from {p = 0,q = 0). Hence we must have a precise relation among the slopes of such two 
curves near the origin, i.e. liuig^Q p{q) / q > limp^Q p/q{p)] but since liuiq^Q p[q) / q = a(3'^ and 
liuip^Q p / q{p) = 1//3^(1 — a), the latter inequality simply leads us to conclude that only trivial 
intersection point are possible for 4/3^a(l — a) < 1. Therefore, we can resume all these results 
in the following 

Proposition 3. In the thermodynamic limit, it exists and it is unique the replica symmetric 
pressure of the bipartite spin glass, given by Ii66\) . Furthermore the region of the {a, /3) plane 
such that the i67\) and /i68\) have only the trivial solutions is characterized by 4/3^a(l — a) < 1. 

Remark 5. In fact for 4/3'^a(l — a)<l the minmax is obtained for q,p = 0, and, as it is easily 
seen, the pressure i66\) reduces to the annealed one ^9^ , that coincides with the true pressure of 
the model in the thermodynamic limit in such a region. 

Furthermore, bearing in mind (j66p . together with (f67j) and (j68p . it is a remarkable result that 
the replica symmetric free energy of the bipartite model is given by the convex combination of 
two monoparty spin glasses, at different temperatures, exactly as happens in the ferromagnetic 
case. This is clarified by the following 

Proposition 4. Choosing f3' = P\/2oi^J~^^ and /3" = f3y^2{l — a)yj'^^, we have 

Ans{a,^) = aAi^iP') + (1 - a)ylgf (/?"), (69) 
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while, with the different scaling of the inverse temperatures /?' = /3\/2ayp and fi" = /3y^2(l — a)y 
we have 

Ansioi, fi) - AA{a, /?) = a (/?') - (/?')) + (!-«) (^|f (/3") - • (70) 

Proof. The proof follows by a straight calculation. If we take for the two monoparty model 
the two different inverse temperatures /3' = /3\/2aa^ and /3" = /3y^2(l — a)a? , with a a free 
parameter to be determined at the end, we have 

ARs{a, f3) = In 2 + oEg log cosh [g|3^J2{l - a)p) + 

+ (1 - a)Eg log cosh(5/3y2^) + a(l - a)f3'^{l - q){l - p) 

= In 2 + oEg log cosh [gl3"a^y/p) + 

+ (1 - a)E<,logcosh(^V^) + Va(l - a)/3'/3"(l - -j)) 

< a ^log 2 + Eg log cosh {g^"a^^/p) + ^a^(l - pf^ 

+ (1-a) flog2 + E„logcoshf5^V^VT^(l-^)'V (71) 



Now it is easily seen that in the last formula we get an equality with = . Following exactly 
the same path of the previous part of the work (dealing with ferromagnetic models) we recover 
(|70p with the choice = |. We stress that this last value of a is more meaningful, in the sense 
that in this case the two monoparty models are trivially independent and separated, with the 
order parameters given by usual self consistency relations for the SK model: 

q = Egtanh^(5r/3'VI), 
p = Egtanh^(5r/3"\/F), 
in perfect analogy with the ferromagnetic case. □ 

Lastly, as it is well known, a theory with no overlap fluctuations allowed may not hold at 
low temperatures and we want to report about its properties in the limit /3 — )• oo to check the 
stability of the replica symmetric ansatz. We will concern about the ground state energy crs 
and its associated entropy sr5, defined by 

eRs{a) = lim 5/3AH5(a,/3) = lim iij5(a,/3)//3, (72) 

/3— >oo 13— ^oo 

SRs{a) = lim. {ARs{a,P) - Pd/^ARsia,^)). (73) 

p—^oo 

First of all, from the self-consistency equations ([671 l68 p . through a long but straightforward calcu- 
lation, we can compute the low temperature limit for the order parameters, q(a, f^),p{a, f^) — )• 1, 
together with the rates they approach to their limit value, /?(1 — g(a, /?)) — )• (7r(l — a))~^^'^ , 
p{l-p{a,P)) (7ra)-V2. Then, bearing in mind the explicit form of the pressure of the model 
in the replica symmetric regime (I66p . we derive the following expressions for the ground state 
energy and the entropy: 



eRsia) = Y^^^, (74) 
SRsia) = (l- Va(l-a)) • (75) 
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Notice that the entropy sptsict) is strictly less than zero for every a G (0,1), that is a typical 
feature of the replica symmetric ansatz for glassy systems. Therefore, the true solution of the 
model must involve replica symmetry breaking. Furthermore, it is a concave function of a, and 
assume its maximum value srs = — 1/vr in a = 1/2, i.e. the balanced bipartite, and its minimum 
value sjis = —2/11 at the ending points a = 0, 1, when the size of one party is negligible in the 
thermodynamic limit. Analogously, for the ground state energy we find in the perfectly balanced 
case ejis{l/2) = 1/2-^7?, and at the extrema of the definition interval of a 6/^5(0) = 6/^5(1) = 0. 

3.3 Constraints 

While the order parameters for simple models (as the bipartite CW) are self-averaging, frus- 
trated systems are expected to show the replica symmetry breaking phenomenon [28||18|. which 
ultimately inhibits such a self-averaging properties for ((712)1(^12)- As a consequence a certain 
interest for the constraints to free overlap fiuctuations raised in the past [17] [16] [T] [5] (and re- 
cently has been deeply connected to ultrametricity [29j [3j) which motivate us to work out the 
same constraints even in bipartite models. 

To fulfil this task the first step is obtaining an explicit expression for the internal energy density 
(which is self- averaging [11]). 

Theorem 5. The following expression for the internal energy density of the bipartite spin glass 
model holds in the thermodynamic limit 



As the proof can be achieved by direct evaluation, we skip it and turn to the constraints: 
Starting with the linear identities we state the following 

Proposition 5. In the thermodynamic limit, and (3 almost-everywhere, the following generaliza- 
tion of the linear overlap constraints holds for the bipartite spin glass 



Proof. Let us address our task by looking at the /3 streaming of the internal energy density, once 
expressed via {qi2Pi2)', hi a nutshell, physically, we obtain these constraints by imposing that 
such a response can not diverge, neither in the thermodynamic limit: 




(76) 



{(I12P12) - 4(^12^12923^23) + 3(gi2Pl2g34P34) = 0. 



(77) 




(78) 



(79) 



now we use Wick theorem on ^ and introducing the overlaps we have 




(80) 



(Pl2Q'l2Pl3gi3) + (Pl2g'l2P34g34) + {pqi2Pl2 " {Pl2qi2Pl3qi3) 

{pi2qi2Pi3qi3) + (^12912^34934) - {pqi2Pi2) + {pi2qi2P3iq3i) 
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The several cancelations leave the following remaining terms 



dl3{qi2Pl2) = N2\{qi2Pl2) - ^{<ll2Pl2q2ZP2z) + Hqi2Pl2q?,iP?,i)j (81) 

and, again in the thermodynamic limit, the thesis is proved. □ 

Proposition 6. In the thermodynamic limit, and in (^-average, the following generalization of 
the quadratic Ghirlanda-Guerra relations holds for the bipartite spin glass 

{qi2Pi2q23P23) = ^(gi2Pi2> + ^{quPuf , (82) 
{qi2Pi2q34P34) = -^(912^12) + '^{qi2Pi2f- (83) 



Proof. The idea is to impose, in the thermodynamic limit, the self- averaging of the internal 
energy (i.e. (e^(a,/3)) — (e(a,/3))^ = 0); we obtain a rest that must be set to zero and gives the 
quadratic control. Starting from 

E(e2(a,/3)) = ^j^^ ^ ^^s ^Y.^^^^^'^^^'^i^jM'^j^'''^^ 

with a calculation perfectly analogous to the one performed in the proof of Proposition [5] and 



comparing with the former relations, we get the linear system 

= {qi2Pl2) + 6(gi2Pl2934P34) " 6(^12^12523^23) " (912^12)^, (84) 

= (gi2Pl2) -4(gi2Pl2923P23) +3(gi2Pl2g34P34), (85) 

whose solutions give exactly the expressions reported in Proposition [6l □ 



4 Conclusion 

In this paper we analyzed the equilibrium behavior of bipartite spin systems (interacting both 
with ferromagnetic or with spin glass couplings) trough statistical mechanics; these systems are 
made of by two different subsets of spins (a priori of different nature |9][15J), for the sake of 
clearness each one interacting with the other, but with no self-interactions. For the former class 
trough several techniques, among which our mechanical analogy of the interpolation method, 
early developed in |17] and successfully investigated in |l][15j[8]. we have seen that the thermo- 
dynamic limit of the pressure does exist and it is unique and we gave its explicit expression in 
a constructive way via a minimax principle. Further, when introducing the Burger's equation 
for the velocity field in our interpretation of the interpolating scheme, our method automati- 
cally "chooses" the correct order parameter, which turns out to be a linear combination of the 
magnetizations of the two subsystems with different signs, so to convert the minimax variational 
principle in a standard extremization procedure. Noticing that the same structure can be re- 
covered for many other models of greater interest, like bipartite spin glasses, we went over and 
analyzed even the latter. 

For these models we have studied both the annealing and the replica symmetric approximation 
(the latter trough the double stochastic stability technique recently developed in ^) which al- 
lowed us to give an explicit expression for the free energy and to discover and discuss the same 
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minimax principle of the bipartite ferromagnets. 

Furthermore, we evaluated the replica symmetric observable in the low temperature limit con- 
firming the classical vision about the need for a broken replica symmetry scheme: one step 
forward in this sense, by studying the properties of the internal energy, we derived all the clas- 
sical constraints to the free overlap fluctuations (suitable obtained for these systems) and we 
worked out a picture of their criticality to conclude the investigation. 

Future works on these subject should be addressed toward a complete full replica symmetry 
broken picture and to a systematic exploration of the multi-partite equilibria. 
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